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Abstract 


It has long been known that many multicomponent alloys 
exhibit surface segregation, in which the composition at or near 
a surface may be substantially different from that of the bulk. 

A number of phenomenological explanations for this tendency 
have been suggested, involving, among other things, differences 
among the components' surface energies, molar volumes, and 
heats of solution. In this work, we investigate via atomistic 
Monte Carlo computer simulation the problem of surface 
segregation for a variety of bcc binary metal alloys, and 
compute segregation profiles (that is, profiles of composition as 
a function of depth) in these materials. Our simulations are 
performed using an atomistic Monte Carlo procedure that 
makes use of the semi-empirical energy method of Bozzolo, 
Ferrante and Smith (BFS) that has been shown to provide 
accurate information on the bulk and surface structure of alloys. 
We formulate a framework for the understanding of this 
phenomenon that makes use of the BFS method's division 
of the energy into "strain" and "chemical" components, and 
compare this picture with available experimental results. 


Introduction 


The composition of metal alloy surfaces is often different from 
that of the bulk. 

Some alloys exhibit surface segregation, where one or more 
species reside preferentially at or near the surface. 

A detailed understanding of this behavior is necessary to 
correctly model such phenomena as adhesion or catalysis. 

Several phenomenological approaches to the problem have been 
put forward, falling into two broad categories: 

Thermodynamic approaches, where the equilibrium 
distribution of chemical species is computed. 

Atomistic approaches, where the tendency of a species to 
segregate is determined by computation of the energies of 
single atoms of that species in bulk and surface 
environments. 


We present results that combine aspects of both approaches: 

Heats of segregation based on total energy calculations 
using the quantum-approximate BFS energy method. 

Equilibrium distribution of species (i.e. segregation 
profiles) computed via atomistic Monte Carlo simulation. 

We attempt to determine whether the heats of segregation 
alone are sufficient to accurately characterize the 
segregation behavior of binary alloys in dilute and non- 
dilute cases. 


Segregation is usually considered to be the result of a 
competition between two effects: 

The effect of size mismatch among the constituent 
chemical species. 

The effect of breaking chemical bonds due to the reduced 
coordination at the surface. 

In this work, we describe the effects of size mismatch in terms 
of the "strain" component of the energy, characterized by the 
degree of mismatch between the lattice constants of the two 
chemical species in a binary alloy. 

We describe the effects of breaking chemical bonds in terms of 
the difference in "chemical" energies of bulk and surface atoms, 
characterized by the chemical identities of these atoms' nearest 
and next-nearest neighbors. 


Method 


The heat of segregation is defined as the difference in energies 
of two otherwise-identical computational cells having a single 
impurity atom in the bulk and in the surface layer. 

Energetics are computed using the BFS quantum-approximate 
energy method. 

The BFS method is especially suitable for analysis of this 
problem, since it explicitly divides the total energy into strain 
and chemical contributions. The strain energy is computed 
under the assumption that each atom's neighbors are in their 
actual positions, but are of the same chemical species as the 
atom itself. The chemical energy is computed under the 
assumption that the neighbor atoms maintain their actual 
chemical identities, but are located at crystallographic positions 
appropriate to a perfect crystal. 

The method has been used previously to compute heats of 
segregation for fee binaries that were in excellent agreement 
with experiment. 


Dilute-limit heats of segregation 

Heats of segregation of single substitutional impurities were 
computed for all binary combinations of Fe,Cr,V,Mo,W,Ta and 
Nb. 

Heats were computed for the 001, Oil and 1 1 1 surfaces, with 
and without atomistic relaxation. 

Segregation trends were characterized as follows: 

E se g = E(surface impurity) - E(bulk impurity), so systems with 
negative heats of segregation will segregate. 

The segregation trend is characterized as strong or weak 
depending on whether the absolute value of the heat of 
segregation is greater than or less than 0.5 eV. 

The trend is characterized as strain-dominated if 
l^strl > | E c h em |, and chemically-dominated if 

l^strl < I E c heml- 

Results for the 001 surface are shown in Table 1. 
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General trends: 


In the following, the size mismatch is characterized by the 
impurity (aj m p) and host (a^ost) lattice constants. 

Binaries having aj m p > aj 10S t tend to segregate, while those 
having aj m p < ahost tend not to. 

Binaries having aj m p > ahost tend to he strain-driven, while 
those having ajmp < ahost tend to he chemically-driven. 

Most exceptions to the above have only a small degree of size 
mismatch, which tends to diminish the importance of strain 
effects and increase the importance of chemical effects. 

Of the binaries having aj m p > aj 10 st that do not segregate, most 
exhibit strong chemically-driven behavior 

Of the binaries having aj m p < ahost> the majority have 
chemically-dominated behavior. 

While relaxation changes the values of the heats of segregation, 
in only three cases does it change the segregation behavior. 


Conclusions: 


In general, if the impurity is larger than the host, strain effects 
will dominate, and segregation will occur. 

In general, if the impurity is smaller than the host, chemical 
effects will dominate, and segregation will not occur. 

Segregation predictions tend to be more accurate for systems 
having large size mismatch. 

Describing segregation behavior as the result of a competition 
between strain and chemical effects is a useful means of 
categorization. 


Extension to Non-Dilute Cases 


Results computed in the dilute limit may not be appropriate for 
non-dilute alloys: 

Non-dilute alloys will have lattice constants different from 
those of the pure-metal hosts, as used in the dilute calculations, 
so the strain contribution to the heat will change with 
composition. 

The neighbors of an "impurity" atom will not necessarily all be 
of the "host" species, so the chemical contribution will be 
composition-dependent as well. 

The chemically-driven formation of large-scale structures in the 
bulk may add further complications. 

We address the first and third of these issues in this work. 

To address the strain-dependence of the heats, we have 
recomputed these heats in the dilute limit, but over a range of 
lattice constants that spans the impurity and host values, 
without atomistic relaxation. 

This approach provides a "strain-dependent" prediction of 
segregation trends for non-dilute binaries which should be more 
accurate than predictions computed in the dilute limit. 


Results: 


Segregation trends (Table 2) are similar to those seen in the 
initial dilute calculations. If the strain-dependent results are 
summarized in terms of the majority segregation prediction for 
each system at all lattice constants, only one system, Fe(Nb), 
exhibits different behavior in the dilute and strain-dependent 
cases. 


Five of the binary systems showed a "crossover" from non- 
segregating to segregating behavior, or vice versa. Of these, 
four systems showed segregation when the "impurity" species 
lattice constant is significantly larger than that of the alloy, but 
not otherwise. 

General trends are consistent with earlier work: 

Binaries having aj m p > a^ost will segregate, while those having 
ajmp < ahost wil1 not 

Of the binaries having aj m p > a^ost that do not segregate, most 
exhibit strong chemically-driven behavior 
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Additional trends: 


Some systems show a reversal of segregation behavior with 
changing lattice constant (Cr(Fe), Fe(Mo), Cr(Mo), Fe(Nb) and 
Cr(Nb)). For the two systems with the largest size mismatch, 
as the lattice constant increases toward the (larger) host value, 
the strain energy becomes less important and the segregation 
behavior is increasingly chemically-dominated. 


Segregation Profiles 


Segregation profiles (plots of composition versus depth) were 
computed via Monte Carlo computer simulation. 

Monte Carlo Procedure 

An initial average composition is chosen, and the chemical 
species distributed at random in a 5200-atom bulk 
computational cell. The compositions are 10/90, 25/75, 50/50, 
75/25 and 90/10. 

The optimum lattice constant of the cell is determined by 
minimizing the total energy with respect to dilation of the cell. 

A surface is created at one face of the cell (001 in this case). 

The chemical species are rearranged probabilistically so as to 
minimize the total energy using the following iterative 
procedure: 

A pair of atoms having different chemical species is chosen. 

The total energy of the two atoms and their environments 
(including nearest- and next-nearest-neighbors of both 
atoms) is computed using the BFS energy method. 

The chemical species of the two atoms are exchanged and 
the total energy recomputed. 


The exchange is accepted or rejected, according to the 
Metropolis criterion: 

If the exchange lowers the energy, it is always accepted; 

If the exchange increases the energy, the exchange is 
accepted with probability exp(-AE/kT) where AE is the 
energy change associated with the exchange, k is the 
Boltzmann constant and T the absolute temperature. 

Issues 

If there is a significant degree of segregation, the bulk 
composition will differ from its initial average value. 

In many cases, the average composition away from the surface 
after minimization is well-defined, and can be taken as the 
effective bulk composition. 

In some cases, e.g. when compositional fluctuations are 
comparable in size to the computational cell, the composition 
away from the surface is not well-defined. 

In addition to the strain/chemical competition, there may be a 
tendency toward bulk order that can mask the tendency toward 
segregation. 

It is expected that the Monte Carlo simulations will be "noisier" 
than the simple dilute-case computations; the MC results will 
exhibit a degree of sensitivity to initial conditions, dependence 
on the particular sequence of random numbers used, etc. 


Results 


Segregation trends from the Monte Carlo simulations are shown 
in Table 3. 

In most cases, the surface segregation trends from the Monte 
Carlo procedure agrees with those from the dilute and strain- 
dependent computations. Most systems show majority 
agreement with the strain-dependent results. 

In 16 of the 42 cases there is disagreement between the MC and 
dilute results for at least one composition. Of these, only four 
systems show disagreement at a majority of compositions. Of 
these, the three have small size mismatches. 
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Key: 

Y=Segregates First Column: Prediction from dilute heats at same lattice constant 

N=Does not segregate Second Column: Prediction from Monte Carlo simulation 


Compositions: 
Host/Impurity = 90/10 
75/25 
50/50 


Specifics: 


Because of the large size of the computational cells, a rich set 
of behaviors can be observed. We consider in detail some 
specific binary systems of interest. 

TaFe (Figs. 1-5) 

The TaFe system has a large degree of size mismatch, and is 
strain-dominated in the Fe(Ta) case and chemically-dominated 
in theTa(Fe) case. This system strongly favors the formation of 
an ordered B 2 structure, either on the surface or in the bulk, 

with a surface layer of Ta. This tendency towards ordering is 
so strong that essentially all of the minority species appears in 
the B 2 region, while the rest of the cell appears as a region of 

majority species pure-metal. 


TaNb (Figs. 6-10) 

The TaNb systems has a very small degree of size mismatch. 

Ta segregates to the surface for all compositions. At small 
concentrations of Ta, the Ta segregates to the surface, but a 
complete layer of Ta does not form, as their is insufficient Ta in 
the computational cell. At intermediate concentrations, the Ta 
segregates to the surface, forming a layer of pure Ta. At larger 
concentrations of Ta, the Ta segregates to a region that includes 
the surface layer and several underlying layers, and Nb tends to 
segregate in the bulk. The Ta25Nb75 cell shows a single layer 
of Ta segregated to the surface, and a damped disordered 
modulated structure having a three-layer periodicity. A 
strongly-ordered B 2 structure is not observed. 


NbW (Figs. 11) 


The Nb75W25 system shows a region of pure Nb that includes 
the surface, and a region of bulk-segregated W. 

VW (Figs. 12-13) 

The VW system shows how competition between bulk and 
surface segregation can lead to varying results. 

V50W50 

Bulk and surface segregation of W. 

V10W90 

Surface segregation of V. 
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Conclusions 


Segregation behavior from dilute heats of segregation, strain- 
dependent heats, and Monte Carlo simulation are in good 
agreement. While there is not a great deal of experimental data 
for the systems studied, the excellent agreement between 
experiment and earlier segregation predictions for fee binaries 
(ref) using the BFS method suggests that this method may 
provide accurate results with bcc systems as well. 


Somewhat surprisingly, even the heat of segregation computed 
in the dilute limit is a good indicator of the segregation 
behavior. 

The good agreement between the predictions of the dilute-limit 
heats and the strain-dependent ones means that little additional 
information is gained from the strain-dependent calculations. 

To obtain a level of detail comparable to that provided by the 
strain-dependent calculations, heats of segregation can be 
calculated using computational cells that have been equilibrated 
using the Monte Carlo scheme described above. 

The MC calculations show a rich spectrum of surface and bulk 
behavior. 


